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Abstract 

One  of  the  most  efficient  approaches  known  for  finding  eun  optimal  tour 
of  the  asymmetric  traveling  salesman  problem  (ATSP)  is  branch-and-bound 
(BnB)  subtour  elimination.  For  two  decades,  expert  opinion  has  been  di¬ 
vided  over  whether  the  expected  complexity  of  the  ATSP  under  BnB  subtour 
elimination  is  polynomial  or  exponential  in  the  number  of  cities.  We  show 
that  the  argument  for  polynomial  expected  complexity  does  not  hold. 


*This  paper  is  to  appear  in  Operations  Research  Letters. 

^The  author  was  partially  supported  by  DARPA  Contract,  #MDA972-94-2-0010. 


1  Introduction 


Given  n  cities  {1, 2,  •  •  • ,  n}  and  a  cost  matrix  (c,j)  which  defines  the  cost  between 
each  pair  of  cities,  the  traveling  salesman  problem  (TSP)  [8]  is  to  find  a  minimum- 
cost  tour  that  visits  each  city  once  and  returns  to  the  starting  city.  When  the 
cost  Cjj  from  city  i  to  city  j  is  not  necessarily  equal  to  that  from  city  j  to  city  i, 
the  problem  is  the  asymmetric  TSP  (ATSP).  It  is  well  known  that  the  ATSP  is 
iVP-hard  [5]. 

One  of  the  most  efficient  approaches  for  optimally  solving  the  ATSP  is  branch- 
and-bound  subtour  elimination  [2,  3]  using  the  assignment  problem  as  a  lower- 
bound  function.  The  solution  of  an  assignment  problem  (AP)  [11]  is  either  a 
complete  tour  or  a  collection  of  disjoint  subtours,  which  is  a  relaxation  of  the 
ATSP  and  can  be  solved  in  O(n^)  time.  Branch-and-bound  (BnB)  [2,  9]  is  one  of 
the  most  general  and  efficient  algorithms  for  finding  the  exact  solutions  of  most 
combinatorial  optimization  problems.  The  algorithm  first  solves  the  AP  for  the  n 
cities.  If  the  solution  is  not  a  tour,  then  the  problem  is  decomposed  or  expanded 
into  subproblems  by  eliminating  one  of  the  subtours  in  the  solution.  A  subproblem 
is  then  selected  and  the  above  process  is  repeated  until  each  subproblem  is  solved, 
i.e.,  its  AP  solution  is  a  tour,  or  until  the  costs  of  the  AP  solutions  of  the  unsolved 
subproblems  are  greater  than  or  equal  to  the  cost  of  the  best  tour  obtained  so  far. 

Bellmore  and  Malone  [3]  argued  that  the  ATSP  can  be  solved  in  polynomial  ex¬ 
pected  time  using  BnB  subtour  elimination.  Their  argument  treated  the  process  of 
BnB  subtour  elimination  as  a  statistical  experiment,  with  the  ^-th  trial  correspond¬ 
ing  to  the  exploration  of  the  z-th  selected  subproblem,  and  success  corresponding 
to  the  z’-th  selected  subproblem  producing  a  tour  of  mininmmn  cost.  Let  the  se¬ 
quence  of  subproblems  selected  by  BnB  subtour  elimination  be  Xq?  A"!,  X2,  A3,  •  •  •, 
and  Pi  be  the  probability  that  the  AP  solution  of  X,  is  a  tour,  for  z  =  0, 1, 2,  •  •  •. 
They  began  with  an  observation  that  the  probability  po  of  the  AP  solution  to 
the  original  problem  being  a  complete  tour  is  approximately  e/n,  where  n  is  the 
number  of  cities.  They  then  assumed  that  the  event  that  the  AP  of  Xi  is  a  tour 
is  independent  of  whether  the  AP  of  Xj  is  a  tour  for  z  ^  j,  and  that  p,  >  po  for 
i  >  1-  Under  these  assumptions,  the  expected  number  of  expanded  subproblems 
is 

00  z— 1  00 

Yj  11(1  “  Pi)  ^  Y  *Po(l  -  =  1/po  =  0{n) 

z=l  z=l 

for  large  n.  Thus  the  expected  running  time  is  O(n^),  as  the  APs  can  be  solved 
in  no  more  than  O(n^)  time  [11].  Their  experiments  show  an  expected 

running  time  for  10  <  n  <  80.  Additional  experiments  reported  by  Smith  et  al.  [14] 
show  that  the  expected  running  time  is  for  30  <  n  <  200.  Smith  [13]  also 

argued,  under  many  assumptions,  that  the  expected  complexity  is  0(n^ln(n)). 


Lenstra  and  Rinnooy  Kan  [10]  pointed  out  that  Bellmore  and  Malone’s  assump¬ 
tions  may  not  be  valid.  In  particular,  two  trials  may  not  be  independent  of  each 
other,  and  pi  may  not  be  greater  than  po,  for  i  >  1.  They  also  pointed  out  that 
to  have  an  expected  complexity  of  0(n‘’),  for  some  constant  c  >  0,  the  following 
two  conditions  must  hold.  First,  the  AP  solution  of  Xq,  the  original  problem, 
must  have  a  probability  po  =  0{n~^)  of  being  a  tour.  Second,  the  number  of  ex¬ 
panded  subproblems  whose  probability  of  yielding  a  tour  is  less  than  po,  must  be 
a  constant.  In  addition,  Balas  et  al.  [1,  2]  carried  out  statistical  analysis  on  the 
experimental  data  obtained  from  three  efficient  implementations  of  BnB  subtour 
elimination.  Their  study  concluded  that  over  the  range  of  40  <  n  <  325,  the 
performance  of  BnB  subtour  elimination  on  randomly  generated  ATSPs  can  be 
described  almost  equally  well  by  a  polynomial  function  (ctn^),  a  superpolynomial 
function  or  an  exponential  function  (ae^”). 

In  short,  whether  the  expected  complexity  of  the  ATSP  under  BnB  subtour 
elimination  is  a  polynomial  or  an  exponential  function  of  the  number  of  cities  is 
an  open  question,  on  which  the  opinion  of  experts  is  divided  [7]. 

In  this  paper,  we  disprove  Bellmore  and  Malone’s  polynomial  argument.  We 
show  that  their  critical  assumption  that  p,  >  po,  for  i  >  1,  does  not  hold  under 
random  ATSPs.  In  particular,  we  show  that,  even  if  the  subproblems  selected 
by  BnB  subtour  elimination  are  assumed  to  be  independent  of  each  other,  the 
algorithm  expands  more  than  ln(n)  number  of  subproblems  Xi  for  which  p,  <  po. 
This  result  supports  and  finalizes  Lenstra  and  Rinnooy  Kan’s  argument. 

In  Section  2,  we  briefly  discuss  the  BnB  subtour  elimination  algorithm  and  the 
analytic  model  we  use.  In  Section  3,  we  disprove  the  polynomial  argument.  Our 
conclusion  appears  in  Section  4. 


2  Preliminaries 

Let  V  be  the  set  of  n  cities,  and  C  =  (c,j)  be  a  cost  matrix  that  specifies  costs 
between  all  city  pairs.  Let  n„  be  the  set  of  all  permutations  of  V,  which  defines 
all  possible  solutions  to  the  assignment  problem  [11]  on  V.  A  permutation  tt  G  n„ 
specifies  an  assignment  of  the  cities,  in  which  city  T(i)  is  assigned  to  city  i,  or  7r{i)  is 
the  successor  of  city  i.  Thus,  under  tt,  a  chain  1  —  7r(l)  —  7r(7r(l)) - 7r---;r(l)  —  1 

forms  a  subtour  if  A:  <  n,  or  a  complete  tour  if  ^  =  n.  When  A;  =  n,  tt  is  called  a 
cyclic  permutation.  Let  11*  C  Iln  be  the  set  of  all  cyclic  permutations  of  V.  For  it  € 
n,  define  /(tt,  C)  =  ^^(i)  to  be  the  cost  of  tt  under  cost  matrix  C.  The  ATSP 

is  to  find  a  cyc/ic permutation  ir*  G  11*,  such  that  /(tt*,  C)  =  min{/(7r,  G  H*}. 
Similarly,  the  assignment  problem  (AP)  seeks  a  permutation  ir'  G  Iln,  such  that 
f{ir',C)  =  min{/(x,  C')|7r  G  n„}.  In  other  words,  the  solution  to  the  AP  is  a 


collection  of  disjoint  subtours.  If  the  AP  solution  happens  to  be  a  complete  tour, 
it  is  the  solution  to  the  ATSP  as  well. 

BnB  subtour  elimination  first  solves  the  AP  for  all  n  cities,  in  time  O(n^)  [11]. 
If  the  solution  is  not  a  complete  tour,  then  the  problem  is  decomposed  into  sub¬ 
problems  by  excluding  some  of  the  arcs  in  a  subtour,  which  eliminates  that  subtour. 
Which  subtour  to  choose,  and  how  to  eliminate  a  subtour,  constitute  the  branching 
rules.  One  heuristic  is  to  select  a  subtour  with  a  minimum  number  of  arcs  that 
axe  not  in  the  included  axe  set  [4],  so  that  the  number  of  subproblems  generated 
from  a  problem  is  minimized,  as  excluding  one  arc  that  is  not  included  is  suffi¬ 
cient  for  breaking  a  subtour.  After  a  subtour  is  chosen,  the  subproblem  should 
be  decomposed  such  that  no  duplicate  subproblems  axe  generated,  minimizing  the 
total  number  of  subproblems  generated.  One  such  branching  rule  was  contributed 
by  Caxpaneto  and  Toth  [4] .  Their  rule  selects  a  subtour  with  the  minimum  num¬ 
ber  of  arcs  that  axe  not  in  the  included  arc  set,  and  decomposes  the  subproblem 
as  follows.  Let  E  be  the  excluded  axe  set,  and  I  be  the  included  arc  set  of  the 
problem  to  be  decomposed.  Assume  that  there  are  t  arcs  of  the  selected  subtour, 
{x\,X2,  -  •  ■  that  are  not  in  I.  The  rule  decomposes  the  problem  into  t  chil¬ 
dren,  with  the  A:-th  one  having  excluded  arc  set  Ek  and  included  arc  set  Ik,  such 
that 


Ek  =  E\J  {x*}, 

h  =/U{®i,-- 


1  U  =  l,2,' 


,t. 


(1) 


Since  Xk  is  an  excluded  axe  of  the  fe-th  subproblem  generated,  Xk  G  Ek,  and  it  is 
an  included  arc  of  the  k  -|-  1-st  subproblem,  Xk  €  Ik+i,  a  tour  obtained  from  the 
k-th.  subproblem  does  not  contain  arc  Xk,  but  a  tour  obtained  from  the  k  +  1-st 
subproblem  must  have  arc  Xk.  Thus  a  tour  from  the  A:-th  subproblem  cannot  be 
generated  from  the  k  -|-  1-st  one,  and  vice  versa.  Therefore,  the  state  space  of  the 
ATSP  under  BnB  subtour  elimination  can  be  represented  by  a  search  tree  without 
duplicates.  This  state-space  tree  can  be  explored  by  best-first  search  or  depth-first 
search. 

Following  previous  research  [6,  7],  in  the  following  analysis  we  use  a  random 
cost  matrix  {cij)  whose  elements  axe  independently  and  uniformly  chosen  from  the 
unit  interval  [0, 1],  and  we  do  not  enforce  the  symmetry  (c,j  =  cj,,)  or  the  triangle 
inequality  {cij+Cj^k  >  Ci.fc)-  Since  we  are  interested  in  finding  a  cyclic  permutation 
of  the  cities,  we  can  set  Ci,i  =  oo  without  loss  of  generality,  for  ^  =  1, 2,  •  •  • ,  n.  The 
AP  of  a  cost  matrix  with  Ci^i  =  oo  is  called  a  modified  AP,  and  a  permutation  with 
Ci^i  =  oo  is  called  a  feasible  permutation  [3]. 


Lemma  2.1  [3]  Given  n  cities,  there  are  asymptotically  \n\je  0.5J  feasible 

permutations  of  the  cities.  □ 
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Lemma  2.2  [3,  7]  Given  an  n  xn  random  matrix,  the  probability  that  a  modified 
AP  solution  is  a  cyclic  permutation  is  asymptotically  e/n.  □ 

Lemma  2.3  Given  an  n  x  n  random  matrix,  the  expected  number  of  subtours  in 
a  modified  AP  solution  is  asymptotically  less  than  ln(n),  and  the  expected  number 
of  arcs  in  a  subtour  of  the  solution  is  greater  than  n/ln(n). 

Proof.  Solving  the  modified  AP  of  a  random  cost  matrix  is  equivalent  to  ran¬ 
domly  selecting  a  feasible  permutation.  This  is  proved  as  follows.  We  first  show 
that  all  feasible  permutations  are  equally  likely  to  have  the  minimum  cost.  Let  n„ 
be  the  set  of  all  feasible  permutations  of  n  cities.  Arbitrarily  partition  n„  into  two 
subsets  and  i.e.,  n„  =  II^  U 11^^)  and  n  =  0.  Then  arbitrarily 
select  one  permutation  tti  from  and  another  one  1^2  from  Without  loss 
of  generality,  assume  that  iti  and  1^2  have  k  arcs  (city  pairs)  not  in  common,  and 
n  -  k  arcs  in  common,  where  A;  =  1, 2,  •  •  • ,  n.  For  a  given  cost  matrix  C,  the  cost 
/(tti)  of  Ti  and  the  cost  /(X2)  of  X2  can  be  written  as 

/(’Ti)  =  +  Y,  =  Y^i+  Y 

*=1  z=l  z=fc+l 

where  e,-  and  s'-  are  different  arc  costs  in  tti  and  7r2.  Obviously,  whether  /(tti)  < 
/(7r2)  or  /(tti)  >  f{ir2)  depends  only  on  the  different  arc  costs,  £,  and  ^i- 
Since  £,•  and  £(•  are  independent  and  identically  distributed  (i.i.d.)  random  vari¬ 
ables,  then  the  probability  of  e,  >  s',  is  the  same  as  the  probability  of 

YlLsi  <  ZLiSi-  In  other  words,  each  of  the  two  permutations  tti  €  11^^^  and 
7r2  €  has  the  same  probability  of  being  smaller.  Consequently,  with  respect  to 
the  permutations  in  11^^),  all  permutations  in  11^^)  are  equally  likely  to  have  smaller 
costs.  Because  the  subsets  11^^^  and  11^^^  were  chosen  arbitrarily,  all  feasible  per¬ 
mutations  are  thus  equally  likely  to  have  the  miniTmim  cost. 

Furthermore,  it  is  known  that  the  expected  number  of  subtours  in  a  permuta¬ 
tion  asymptotically  approaches  ln(n)  [12].  Thus  the  expected  number  of  arcs  in  a 
subtour  of  an  AP  solution  is  0(n/ ln{n)).  □ 

3  On  the  Polynomial  Argument 

Let  us  assume  the  validity  of  Bellmore  and  Malone’s  assumption  that  two  sub¬ 
problems  chosen  by  BnB  subtour  elimination  are  independent  of  each  other.  Then 
their  polynomial  argument  critically  depends  on  the  assumption  that  a  constant 
number  of  subproblems  Xi  are  expanded  whose  probabilities  pi  are  less  than  po- 
In  this  section,  we  show  that  this  argument  does  not  hold  even  if  the  independence 
assumption  is  granted.  We  first  prove  the  following  lemma. 
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(a)  q(n9s4)  when  there  is  no  included  arc. 


(b)  q(n^,t)  when  there  are  included  arcs. 


Figure  1:  Probability  q(n,s,t)  that  the  AP  solution  is  a  tour 


Lemma  3.1  Given  an  n  x  n  random  cost  matrix,  let  q(n,s,t)  be  the  probability 
that  the  solution  of  a  modified  AP,  which  has  s  <  n  excluded  arcs  and  t  included 
arcs,  is  a  tour.  Then  q(n,s,t)  is  asymptotically 

e/n  —  o{l/n)  <  q{n,s,t)  <  efn  +  o(lfn),  when  t  =  0  (2) 

9(n,s,f)  =  +  o(l/n),  whent>0  (3) 

where  \,  1  <  X  <  e,  is  a  constant. 

Proof.  See  Appendix. 

Since  Lemma  3.1  is  important  to  our  analysis,  we  verified  its  correctness  by  ex¬ 
perimentally  solving  10, 000  randomly  generated  APs  for  each  n  €  {10, 20, 30,  •  •  • ,  100}. 
Our  results  are  presented  in  Figure  1.  Figure  1(a)  shows  that  the  experiments  sup¬ 
port  (2)  when  s  =  n/2  excluded  arcs  and  t  =  0  included  arcs  are  used.  Figure  1(b) 
shows  that  the  experiments  support  (3)  when  s  =  n/5  excluded  arcs  and  t  =  nf2 
included  arcs  are  employed. 

A  subproblem  Xi,  which  is  selected  by  BnB  subtour  elimination  and  which 
generates  an  optimal  tour,  must  be  a  minimum-cost  leaf  node  in  the  search  tree. 
Thus,  the  probability  of  X,  generating  an  optimal  tour  cannot  be  greater  than  the 
probability  of  Xi  being  a  leaf  node.  Note  that  Xj  can  generate  a  complete  tour  and 
become  a  leaf  node  only  when  all  its  ancestors  in  the  search  tree  do  not  generate 
complete  tours  but  X,  does. 

Lemma  3.2  Suppose  that  two  nodes  chosen  by  BnB  subtour  elimination  are  in¬ 
dependent  of  each  other.  Let  p  be  the  probability  that  a  non-root  node  in  the  search 
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tree  is  a  leaf,  and  in  'particular  let  pq  be  the  probability  that  the  root  is  a  leaf.  For 
a  non-root  node  "with  t  included  arcs,  there  exists  a  constant  Q  <  8  <1  —  Ife  such 
that,  if  t  <  8n,  then  p  <  po,  where  n  is  the  number  of  cities. 


Proof.  Consider  a  node  K  of  a  search  tree  which  is  generated  by  equation 
(1).  Let  the  number  of  included  arcs  of  F  be  t  and  the  number  of  excluded  axes 
of  Y  be  5,  which  is  also  the  depth  of  the  node  in  the  search  tree.  Denote  the 
nodes  on  the  path  from  the  root  to  the  node  Y  as  where 

Yo  is  the  root.  From  equation  (1),  Yi  has  i  excluded  arcs.  In  addition,  let  Yi 
have  ti  included  arcs.  By  equation  (1),  0  =  to  ■<  ti  <  t2  <•••  <  ts-i  <  t. 
It  is  only  when  none  of  io, Fi,  •  •  • , Fs_i  is  a  leaf  node  that  F  exists  in  the  tree. 
Therefore,  the  probability  q(n,s,t)  of  F’s  AP  solution  yielding  a  complete  tour  is 
equal  to  the  probability  that  all  its  ancestors  do  not  generate  complete  tours  and 
F  produces  a  complete  tour.  The  probability  that  F’s  parent  F-i  does  not  yield 
a  complete  tour  is  (1  —  q{n,  s  —  1,  and  that  F’s  grandparent  F-2  does  not  is 

(1  —  q(n,  s  —  2,ts-2)),  and  so  on.  Consequently,  by  the  independence  assumption, 
the  probability  p  that  F  exists  and  is  a  leaf  is  then 


S—1 

p  =  9(n,s,t)n(l  -q{n,i,ti)). 

i=Q 

By  Lemma  3.1,  we  have 


"  =  (;r^  +  “(V"))  n  (i  - 

where  Aq  >  Ai  >  •  •  •  >  A^.i  >  A 
that 


n-ti^ 

>  1  are  constants.  It  can  be  shown  by  induction 


^  A,-  _  A^s 

^n-ti  n-f' 


1 

where  A'  =  -  ^  A,- 

2  i=0 


and 


t'  = 


^iti 

EI=o  ■ 


Obviously, 


0<t'  <t.  (4) 

The  probability  p  can  be  further  written  as 

p  =  (5) 

We  now  show  the  lemma  by  contradiction.  Assume  that  the  lemma  does  not  hold, 
namely  p  >  po,  where  po  =  e/n  by  Lemma  2.2.  Let  8  =  (e  —  A)/e.  Since  1  <  A  <  e. 


we  know  that  0  <  ^  <  1  —  1/e.  Then  by  ignoring  the  o(l/n)  term  and  some 
algebra,  it  can  be  shown  from  (5)  that  when  t  <  Sn,  p  >  po  is  equivalent  to 


t'  ">■  n  -\- 


XX'sn 

(e  —  A)ra  —  et 


>  n. 


Since  n  >t,we  have  t'  >  t,  which  is  in  contradiction  with  (4).  □ 

In  a  search  tree  with  nodes  generated  according  to  equation  (1),  only  one 
node  on  the  first  depth  of  the  tree  has  no  included  arcs.  When  a  subtour  with  the 
minimum  number  of  arcs  is  chosen,  the  number  of  children  generated  is  0(n/  ln(n)) 
on  average  by  Lemma  2.3.  Thus,  the  nodes  at  the  first  depth  have  t  =  0{n/\n{n)) 
included  arcs.  Asymptotically  t  <  6n,  where  0<^<1  —  1/eisa  constant.  Thus, 
all  nodes  except  one  without  included  arcs  on  the  first  depth  asymptotically  satisfy 
t  <  Sn.  Similarly,  all  nodes  except  one  on  the  i-th  depth  satisfy  t  <  Sn  when  i  is 
no  bigger  than  0(ln(n)).  Now  consider  a  node  that  has  no  included  arcs.  Notice 
that  none  of  the  ancestors  of  the  node  have  included  arcs  either.  By  Lemma  3.1, 
the  probability  of  the  AP  solution  of  the  node  or  one  of  its  ancestors  being  a  tour 
can  be  asymptotically  approximated  as  e/n  by  ignoring  the  o{l/n)  term.  Thus, 
the  probability  p  that  a  node  with  no  included  arcs  exists  and  is  a  leaf  node  on  the 
d-th  depth  of  the  search  tree  can  be  approximated  by  p  =  (e/n)(l  —  e/nY,  which 
is  less  than  e/ n,  the  probability  po  that  the  root  node  is  a  leaf.  Overall,  for  nodes 
at  depth  i  <  no  bigger  than  0(ln(n)),  the  probabilities  that  they  are  leaf  nodes 
are  asymptotically  less  than  the  probability  of  the  root  being  a  leaf. 

Recall  that  the  probability  that  the  AP  solution  of  a  subproblem  selected  by 
BnB  subtour  elimination  will  yield  an  optimal  tour  is  less  than  the  probability 
that  the  node  in  the  search  tree  will  become  a  leaf  node.  If  the  depth  of  the  node 
that  generates  the  optimal  tour  is  greater  than  ln(n),  or  a  node  with  depth  greater 
than  ln(n)  is  expanded,  then  in(n)  nodes  must  be  expanded  whose  probabilities  of 
generating  an  optimal  tour  are  less  than  po.  Consider  the  case  when  the  depth  of 
the  node  generating  an  optimal  tour  is  less  than  ln(n),  and  no  node  at  depth  greater 
than  ln(n)  is  generated.  Assume  further  that  Bellmore  and  Malone’s  argument 
holds,  i.e.,  only  a  polynomial  number  of  nodes  need  to  be  expanded,  which  is 
greater  than  ln(n).  Then  the  probabilities  that  all  these  expanded  nodes  generate 
the  optimal  tour  are  asymptotically  less  than  po,  which  contradicts  the  assumption 
used  in  the  polynomial  argument.  Therefore,  we  have  the  following  result,  which 
disproves  Bellmore  and  Malone’s  polynomial  argument. 


Theorem  3.1  Let  Xq,  X.\.,  X2,  •  •  •  he  the  subproblems  expanded  by  BnB  subtour 
elimination  on  a  random  ATSP,  and  pi  he  the  probability  that  the  modified  AP 
solution  of  Xi  is  a  complete  tour,  for  i  =  0, 1,2,  •••.  BnB  subtour  elimination 
expands  more  than  ln(n)  number  of  subproblems  X,  with  probabilities  pi  <  po  for 
i>l.  U 
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4  Conclusion 


For  two  decades,  the  question  of  whether  the  expected  complexity  of  the  ATSP 
under  BnB  subtour  elimination  is  polynomial  or  exponential  has  remained  open. 
We  have  shown  in  this  paper  that  the  polynomial  argument  is  not  valid. 

The  polynomial  argument  critically  depends  on  the  assumptions  that  for  a 
random  ATSP,  the  subproblems  considered  by  BnB  subtour  elimination  are  inde¬ 
pendent  of  each  other,  and  pi  >  pa,  for  *  >  1,  where  pi  is  the  probability  that  the 
solution  of  the  assignment  problem  to  the  i-th  subproblem  selected  by  the  algo¬ 
rithm  is  a  complete  tour.  We  have  proved  that,  even  if  the  subproblems  selected 
by  BnB  subtour  elimination  are  assumed  to  be  independent  of  each  other,  the 
algorithm  expands  more  than  ln(n)  number  of  subproblems  Xi  for  which  pi  <  po. 


Appendix 

Lemma  3.1  Given  an  n  x  n  random  cost  matrix,  let  q(n,s,t)  be  the  probability 
that  the  solution  of  a  modified  AP,  which  has  s  <  n  excluded  arcs  and  t  included 
arcs,  is  a  tour.  Then  q{n,s,t)  is  asymptotically 

e/n  —  o{l/n)  <  q{n,s,t)  <  e/n  +  o{l/n),  when  t  =  0  (6) 

q(n,s,t)  =  — - |-o(l/n),  when  t  >  0  (7) 

n  —  t 

where  \,  1  <  X  <  e,  is  a  constant. 


Proof.  Denote  by  E  and  I  the  excluded  and  included  sets.  Following  the  same 
argument  used  in  the  proof  of  Lemma  2.3,  it  can  be  easily  shown  that  solving 
the  AP  with  constraints  E  and  I  is  equivalent  to  arbitrarily  selecting  a  feasible 
permutation  from  all  feasible  ones  among  which  there  are  some  cyclic  permuta¬ 
tions.  In  other  words,  the  probability  q(n,  s,  t)  of  the  AP  solution  being  a  cyclic 
permutation  is  equal  to  the  ratio  of  the  total  number  of  cyclic  permutations  to  the 
total  number  of  feasible  permutations  under  the  constraints.  Let  R{n,  s,  t)  be  the 
number  of  cyclic  permutations  and  Q{n,s,t)  the  number  of  feasible  permutations 
with  constraints  E  and  /,  then 


q{n,s,t) 


R{n,s,t) 

Q{n,s,ty 


(8) 


(a)  First  consider  ^(re,  0,  t),  the  case  when  there  exist  some  included  arcs,  but 
no  excluded  arcs.  Assume  initially  that  no  two  arcs  of  I  share  a  common  vertex. 
Consequently,  the  number  of  included  arcs  must  be  less  than  half  of  n,  the  number 
of  cities,  i.e.,  t  <  [n/2J.  For  an  included  arc  (?,  j),  since  i  must  be  assigned  to  j, 


all  out-going  arcs  from  i  and  in-coming  arcs  to  j  except  {i^j)  can  be  ignored.  Thus 
the  AP  with  {i,j)  included  can  be  solved  by  treating  cities  i  and  j  as  one  vertex.  If 
all  arcs  in  I  are  thus  considered  as  cities,  by  Lemma  2.1  there  are  (n  —  t)\le  feasible 
permutations  that  have  no  self  loops,  which  also  exclude  the  feasible  assignments 
of  j  to  i  for  included  arcs  {i,j)  €  I.  We  count  the  number  of  permutations  with 
this  type  of  legal  ‘self  loop’,  namely  i  to  j  and  j  to  i  for  all  {i,j)  6  /,  as  follows. 
For  {i,j)  €  I,  assign  j  to  i,  and  considering  the  other  t  —  1  included  arcs  as  cities, 
there  are  (n  —  2  —  (<  —  l))!/e  permutations,  and  the  total  number  of  choices  for 

{i,j)  is  ^  ^  j-  Similarly,  for  two  arcs  (ij)  and  {i'J'),  assign  j  and  j'  to  i  and 

i'  respectively.  There  are  (n  -  4  -  (i  -  2))!/e  feasible  solutions,  and  the  number 

of  combinations  of  these  two  arcs  is  ^  2 

permutations  with  included  set  I  is 


Thus,  the  total  number  of  feasible 


(5(n,0,t)  =  g  ^  ) 


Obviously, 


Thus,  by  substituting  (11)  into  (9)  we  have 
Q{n,0,t)  <  — 


k=zO 


1 

1 

(9) 

t)\/e. 

it 

(10) 

(n-t)! 
jfc!  ■ 

(11) 

A:=0 

(12) 

Now  consider  the  case  when  there  are  r  cities  which  are  shared  by  two  arcs  in 
I.  Two  arcs  sharing  a  common  vertex  can  be  treated  as  one  arc  connecting  three 
cities,  and  the  shared  vertex  can  be  simply  ignored.  The  problem  is  then  equivalent 
to  one  with  n  —  r  cities  and  t  —  r  included  arcs.  Following  the  previous  reasoning, 
we  have 


(5(n,0,t) 


^  (  t  —  r  \  {n  —  t 

k[  k  )  ^ 


e 


(13) 


Obviously,  inequality  (10)  is  still  valid  in  this  case,  and  the  proof  of  (12)  follows. 
When  there  are  t  included  arcs  and  r  cities  commonly  shared  by  the  included  arcs. 


the  number  of  distinct  cities  supporting  the  included  arcs  is  2t  —  r  which  must  not 
be  larger  than  n,  i.e., 


2t  —  r  <n. 


With  inequality  (14),  it  can  be  shown  that 


Then  the  upper  bound  in  (12)  follows  by  substituting  (15)  into  (9).  Overall,  for  a 
given  included  set  I  with  t  arcs,  we  have 

(n  —  i)!/e  <  (5(n,  0,t)  <  (n  —  <)!.  (16) 

The  number  of  cyclic  permutations  is  simply 

i?(n,  0,t)  =  (n  —  t  —  1)!.  (17) 

Finally,  by  (8),  (16)  and  (17),  we  write 

<  q{n,  0,  t)  <  — (18) 
n  —t  n  —t 

(b)  Now  consider  g(n,s,0),  the  case  when  there  are  some  excluded  arcs,  and 
there  is  no  included  arc.  For  this  case,  we  have 

(5(n,s,0)  =  n!/e  -  (3(n,0,s).  (19) 

This  is  the  difference  between  the  number  of  feasible  solutions  without  any  con¬ 
straints,  which  is  n!/e,  and  the  number  of  feasible  solutions  if  we  include  the 
excluded  arcs  in  the  solutions.  Combining  (19)  with  (16),  we  obtain 

n!/e  —  (n  —  s)!  <  (5(n,s,0)  <  n!/e  —  (n  —  s)!/e.  (20) 

To  compute  i?(n,  s,  0),  we  first  assume  that  no  two  arcs  in  E  have  a  common  ver¬ 
tex.  i?(n,  s,  0)  is  the  number  of  cyclic  permutations  under  the  constraints  of  the 
excluded  set  E.  There  axe  R{n,  0, 0)  cyclic  permutations  without  the  constraints 
of  E.  In  other  words,  we  over-counted  the  number  of  feasible  cyclic  permutations 
by  including  those  arcs  that  were  excluded  before.  The  number  of  cyclic  permu¬ 
tations  with  one  particular  excluded  arc  included  is  R{n,  0, 1),  and  the  number  of 


choices  for  this  arc  is 


.  Similarly,  the  number  of  cyclic  permutations  with  two 
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particular  excluded  arcs  included  is  i?(n,  0,2),  and  the  number  of  combinations  of 
these  two  arcs  is  f  f  ')  •  By  the  principle  of  inclusion  and  exclusion  [12],  we  have 


Obviously,  i?(n,s,0)  is  bounded  by 

R{n,  0, 0)  —  sR{n,  0, 1)  <  R{n,  s,  0)  <  i?(n,  0, 0).  (21) 

Now  we  show  that  (21)  still  holds  when  E  contains  r  cities  shared  by  two  arcs.  The 
upper  bound  holds  since  R{n,  0, 0)  is  the  number  of  cyclic  permutations  without 
any  constraints  posed  by  E.  The  lower  bound  is  valid  because,  by  the  same  reason 
as  used  in  (13),  we  have 


when  there  are  r  common  cities  shared  by  arcs  of  E.  Furthermore,  when  s  <  n. 


Rin,0,k)  <  ^  j  R{n,0,k)  <  ^  j  B(n,0,l). 

e(n  —  1)! 


Using  (17),  (20)  and  (21),  we  write 
e  ((n  —  1)!  —  s(n  —  2)!) 


<  9(n,s,0)  < 


n 


!  —  (n  —  s)!’ 


n\  —  e(n  —  s)! 

By  some  algebra,  we  obtain  the  following, 
e  es  /■  ^ 

n  n(n  _  ^  ^  n  n{n{n  —  1)  •  •  •  (n  —  s  +  1)  —  1) 

When  s  <  n,  we  can  write 


es 


e/n  -  o(l/n)  <  9(n,s,0)  <  e/n  +  o(l/n),  (22) 

which  is  (6). 

(c)  By  Lemma  2.1,  g(n,  0,0)  =  e/n.  Thus  (22)  means  that  asymptotically,  the 

arcs  of  E  have  only  secondary  influence  on  the  probability  that  a  modified  AP 

solution  is  a  tour  when  s  <  n.  Combining  (18)  and  (22),  and  the  fact  that  the 

excluded  and  included  sets  are  disjoint,  we  obtain 

1  c 

<q{n,s,t)<-^  +  o{l/n).  (23) 

Therefore,  there  must  exist  a  constant  A  €  (l,e)  such  that 

^(n,s,t)  = +  o(l/n),  (24) 

which  is  (7).  □ 
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